from pymol import *
cmd.load("Psi_a_126_126-A.cube")
cmd.load("Psi_b_126_126-A.cube")
cmd.load("Psi_a_127_127-A.cube")
cmd.load("Psi_b_127_127-A.cube")
cmd.load("target.mol")
cmd.isomesh("HOMO_A", "Psi_a_126_126-A", -0.02)
cmd.isomesh("HOMO_B", "Psi_b_126_126-A", 0.02)
cmd.isomesh("LUMO_A", "Psi_a_127_127-A", 0.02)
cmd.isomesh("LUMO_B", "Psi_b_127_127-A", -0.02)
cmd.color("blue", "HOMO_A")
cmd.color("red", "HOMO_B")
cmd.color("blue", "LUMO_A")
cmd.color("red", "LUMO_B")
cmd.disable("LUMO_A")
cmd.disable("LUMO_B")
